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Abstract: We prove that SU(N) bosonic Yang-Mills matrix integrals are conver- 
gent for dimension (number of matrices) D > D c . It is already known that D c = 5 
for iV = 2; we prove that D c = 4 for N = 3 and that D c = 3 for N > A. These 
results are consistent with the numerical evaluations of the integrals by Krauth and 
Staudacher. 



Keywords: 



Matrix Models, M(atrix) Theories, INonperturbative Effects. 



Contents 



B 


Introduction 


s 


II- 


Preliminaries 


I 


g. 


£17(2) 


I 


RD 


S77(3) 




& 


iV > 3 




1. 


Discussion 


10 




|6.1| Correlation Functions 


1C 




|6.2| Supersymmetric Integrals 


11 




|6T3| Other Gauge Groups 


11 


0. 


Divergent Matrix Integrals 


12 



1. Introduction 

The discovery of D-branes and the realization of their importance in string theory and 
M-theory has led to a number of exciting conjectures relating these very complicated 
theories to (at least technically) much simpler M(atrix) theories |lj and the IKKT 
model for the type IIB superstring (for reviews see [§,£§]). This has generated re- 
newed interest in supersymmetric Yang-Mills quantum mechanics |], 0, ^], |S| which 
is obtained by the the dimensional reduction of D = 10 supersymmetric Yang-Mills 
gauge theory (SSYM) to 1 remaining space-time dimension. Further dimensional 
reduction to dimensions leads to the supersymmetric Yang-Mills matrix integrals 
which are an important component in calculating the Witten index for the quantum 
mechanics. Both the quantum mechanics and the matrix integrals exist for any di- 
mension D in which the original SSYM exists, ie D = 3, 4, 6, and 10. Dropping the 
supersymmetry requirement leads to the bosonic Yang-Mills matrix integrals which 
exist for any D and are the main subject of this paper. 



A couple of years ago Moore, Nekrasov and Shatashvili [10| found a way of cal 



culating the partition function by deforming the Yang-Mills matrix integrals to a 
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cohomological theory. Their method does not allow the calculation of arbitrary cor- 
relation functions in the original Yang-Mills picture where it is not known how to 
do exact calculations unless the gauge group is SU(2) [JTTJ] . However it is possible 
to do numerical calculations provided the gauge group is not too big. The numeri- 



cally calculated partition functions for small N = 3, 4, and 5 [12, 1 3 1 agree with the 



cohomological calculations in ||10|| . Some correlation functions and eigenvalue distri- 



butions have also been obtained [^3], [14] ; in the case of D = 4 where the fermions 



can be integrated numerically by Monte Carlo many correlation functions have been 



found for values of N up to 48 ||15|| . In the course of their work Krauth and Stau- 
dacher Jl3| also investigated the properties of purely bosonic Yang-Mills integrals 
(these are defined below). It had been believed that the flat directions in the action 
would cause these integrals to diverge (their supersymmetric cousins being saved by 
the Pfaffian arising from the integration of the fermions which vanishes along the 
flat directions). Simple analytic calculations of the partition functions in the case 
of SU(2), and delicate numerical computations for £{7(3), SU(4:) and SU(5) (and 
subsequently for N up to 256 [16[| and other gauge groups as well fi~7[0 showed that 
this is not necessarily the case. Unfortunately up to now an analytic demonstration 
of the convergence of these integrals for SU (N > 2) has been lacking. Our purpose 
here is to provide such a demonstration. 

The bosonic Yang-Mills partition function for gauge group G in D "space-time" 
dimensions 1 dimensionally reduced to is given by 

Zd,g = II / dX ^ E Tr l X ^ X »} 2 C 1 - 1 ) 

where the matrices {X^, /i = 1, . . . D}, which are traceless and hermitian, take values 
in the Lie algebra of G and can be written 

9 
a=l 

The {t a , a = 1, . . . g} are the generators in the fundamental representation satisfying 

Trt a t b = 25 ab (1.3) 
and we shall use / to denote the rank of the Lie algebra. In this paper we will restrict 



ourselves to the groups SU(N). The measure and the integrand in [TTT] are then 
invariant under the SU (N) gauge symmetry 

X„ -> rfXpU, U G SU(N) (1.4) 

1 In the bosonic case there is no requirement of supersymmetry to restrict D so it is possible to 
consider D as a continuous variable by analytic continuation; however in this paper it is to be taken 
strictly as an integer. 
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and the SO(D) symmetry 

QeSO(D). (1.5) 

We will prove that Zd,su(n) is convergent for dimension D > D c and divergent 
for D < D c . It is already known from exact calculation that D c = 5 for N = 2 
(although we will show that our methods reproduce this almost trivially). We prove 
that D c = 4 for N = 3 and that D c = 3 for N > 4. 

The body of this paper is concerned with establishing which integrals converge. 
In section 2 we set up our procedure and establish some results that are useful in 
every case. Section 3 deals with SU(2), section 4 with £77(3), and section 5 with 
SU(N > 3). In Appendix A we show which integrals diverge. We conclude with 
a brief discussion of the implications of our results for the supersymmetric theories 
and other gauge groups. 



2. Preliminaries 



The dangerous regions which might cause the integral |1.1| to diverge are where all 
the commutators almost vanish but the magnitude of goes to infinity. Hence we 
let 

A M = Rx^, Trx^x^ = 1 (2.1) 

where, as from now on, we use the summation convention for Greek indices. Then 
we have 



where 



and 



poo 

Zd,g= / R D9 - x X D , G {R)dR (2.2) 
J o 

D 

Xd,g(R) =]\ dx v 5(1 - Tra^) exp (-R 4 S) (2.3) 



2 

We note that for any finite R the integral Xd,g{R) is bounded by a constant (since 
every term in the argument of the exponential is negative semi-definite) and therefore 
if for large R 

const 

X D , G (R) < -^ r , where a > Dg, (2.5) 

then the partition function Zd,g is finite. Our tactic for proving convergence of Z£,,g 
is therefore to find a bound of the form |2.5| on Xd,g{R)- 
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Now we split the integration region in |2.3| into two 

Tl x : S < (R-^f 



K 2 : S> (R 



-(2-r,)\2 



(2.6) 



where rj is small but positive. We see immediately that the contribution to Id,g{R) 
from TZ 2 is bounded by A\ exp(— R 2r] ) (we will use the capital letters A, B and 
C to denote constants throughout this paper) and thus automatically satisfies 23 



Thus we can confine our efforts to the contribution from TZ\ in which we replace the 
exponential function by unity to get the bound 



X D , G (R) < A\ exp(— i? 2,? ) + X DG [R) 



where 



D 

Since S is a sum of squares it follows that the region TZ[ defined by 



(2.7) 



(2- 



(2.9) 



is larger than the region 1Z± which we can therefore replace in |2.7| by TZ[. 

Now we utilise the SU(N) symmetry to diagonalise x% which we may therefore 
write as x\ = diag(Ai, . . . Ajv). The constraint |2.9| then becomes 

\(.Xi - \j)(x u )ij\ <R- {2 - v) 



(2.10) 



This immediately leads to the generic case which is when the eigenvalues of Xi are 
not degenerate 



I A,- — A,! > e 



(2.11) 



where e is a constant which we may choose but will always be finite. Then all the 
off-diagonal elements of {x u , v — 2, . . . D} are bounded by [2.10| and I D G (R) (|2~ 

D,G\ 

l 9 D en G {R) < A 2 R- {2 -^ D - 1){9 ~ l) (2.12) 



has a contribution l^^R) which is bounded 



where the constant A 2 comes from the integral over all the diagonal elements (and 
does of course depend on e). Note that once the off-diagonal elements are bounded 
by [2.10| and 2.11 all commutators are constrained to be 0(R~^ 2 ~^) although the 
coefficient may be more than 1. Enforcing |2J^ implies constraints on the diagonal 



elements which lower A 2 in |2.12| but does not affect the power of R. 

We note that the case when all eigenvalues are widely separated corresponds 
to the "perturbation" expansion of |Tj|; using |2.5| and |2.12j leads to the criterion 
N > D/(D — 2) for convergence that is given in |16 |. It is fortuitous that this gives 



the correct conditions for convergence of Zb,g because as we shall see X D G is not 
usually the most divergent contribution. 
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3. 577(2) 



In the case of SU{2) each matrix x^ has eigenvalues A M and — A M . However the 
constraint in |2.1| implies that 

2A M A M = 1 (3.1) 

and it follows that there must always be one matrix with eigenvalues of magnitude at 
least {2D)' 2. We can choose this matrix to be x\ and hence by taking e < 2(2D) 2 
we always have the generic case. Hence 

X D ,su(2)(R) < A, exp(- J R 2 ") + A 2 R-* 2 -^ D -V (3.2) 

and it follows immediately that Zd,su(2) is finite for D > 5. We show in Appendix 
A that Z D ^ SU (2) is divergent for smaller D. Of course these results are well known 
because Z D ^u{2) can be calculated exactly. 



4. 577(3) 



From the constraint in 2A it follows that there must always be one matrix with an 
eigenvalue of magnitude at least (3_D)~2; as before we choose this matrix to be x\. 
When all the eigenvalues of X\ are separated by at least e (|2.11|) we get the generic 



contribution [2.12| to the integral. However now we have the new possibility that two 
of the eigenvalues of X\ become degenerate (it is not possible for all three eigenvalues 
to become degenerate provided we choose e < \{?>D)~^) in which case the condition 



2.11| does not apply and we have to proceed more carefully; Id,su(3)(R) is made up 

5(7(3)1 



of a piece where [2.11| applies to all eigenvalues plus a new piece ^jj 9 su(3) (-^) fr° m the 



region of integration where |2.11| is not satisfied by all eigenvalues, 



1, 



First we write 



D,SU(3) 



x x 



(R) 



l^l TT ,JR)+l de9 



-D,su{3y 



D,SU(3) 



(R). 




(4.1) 



(4.2) 



where X\ = diag (£, — £); we know that |2pi/\/3| > (3-D) _ 2 and, because we are 
just interested in the case of degenerate eigenvalues, 2|£| < e. Because only two of 
the three eigenvalues of x\ are degenerate the constraint |2.10| bounds some of the 
off-diagonal elements of £2,...d so we can write these matrices in the form 



x,, 









0{R 



-{2-r,Y 



0(R-^) 




(4.3) 
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where x u is a 2 x 2 traceless hermitian matrix (ie it lives in the su(2) sub-algebra of 
the su(3) algebra inhabited by x v ), and by 0(R^ ( - 2 ^) we mean that the elements 
are bounded by \2.10[ At this stage the off-diagonal elements in the third row and 
column of the x 2 ,...d are innocuous and can be integrated out to get 

z£i m (R) < Bi irw< D -»-*> (f f[ dp,) (7 den 



ttsu(2)(V3pi - e) 2 (V3 Pl + 2 #(1 - (Trx a x a + 2p aPa )) 
#((Tr5 Q x a + 2 PqPq ) - (l - 4 (eR^)~ 



(4.4) 



where 9 denotes the step function, we have included the Vandermonde determinant 
for xi, and flsu(2) is the volume of 577(2). The region 1Z is defined by |£| < e/2 and 



[x,, x v \ 



(4.5) 



Within the region of integration we can bound the factors (3pi ± e) 2 by a constant 
and then, because we are looking for an upper bound, we can drop the constraints 
1^1 < e/2 and |pi| > (4D) _ 2. Doing the integrals we are left with 



l d » 9 S u {3) (R)<B2R- {2 -^ D - 1) Fn 



5(7(2) 



(4.6) 



where 



= (2 - rfiR-^) f du u i-v ( [ Y[dx u ) 5 



R. 



Tr x^x, (4.7) 



Making the rescaling x, = y^u/R} 1 11 ' 2 we find that 

sum = R- {2 ^ )3D/2 f R u^/ 2 )-^ su{2) (u)du. 

Jo 



(4.{ 



As R — > oo we use the results from section 3 to bound the remaining integral giving 



D, 5(7(2) 



< B 3 R-^ 3 °/ 2 , D > 5 

B 4 R-^ 3 °/ 2 log R, D = 4 
BsJr <?-v)QD-i)/2 D = 3 



and so we find that 

X D g su(3)( R ) < B s 



R 2- V 



4(£>-l)+(3£>-5 Dl3 )/2 



(\ogRy D ' 4 



(4.9) 



(4.10) 
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From 12.121 we know that 



(4.11) 



Applying the criterion ^]5| we see that both I^suca) (^) anc ^ ^dsu(3) (^) ma ke a finite 
contribution to Zd,su(3) on ly f° r -D > 4. It is straightforward to check that -Z_d=3 i 5[/(3) 
diverges (see appendix). Thus we conclude that Zd,su(3) is finite for D > 4 only. 

5. JV > 3 

The argument for higher N has the same structure as for SU(3). Again there is a 
generic contribution which is bounded as shown in |2.12| and a degenerate contribu- 
tion. As usual xi can be chosen so that at least one of its eigenvalues has magnitude 
(ND)~2 or more. The difference is that x± can have not only two degenerate but 
2, 3, ... AT — 1 and two or more sets of them. To deal with this we note that whenever 
x\ has some configuration of exactly degenerate eigenvalues there is a sub-algebra of 
su(N) which commutes with X\. (As we saw in section 4 for the case of SU(3) the 
only possibility is to have two degenerate eigenvalues and the sub-algebra is su(2).) 
Suppose that x\ has K sets of degenerate eigenvalues with degeneracies {N ly . . . N K }; 
then the sub-algebra which commutes with x\ is 



K 



H 



su(N k ) 



(5.1) 



k=l 



into 



Now let t a denote those generators of G that do not lie in H, and decompose x^ 

x x = diag { P \I N \ . . . p^I NK } a\ . . . of) + diag (x\, . . . xf , . . . 0) 
x^t = diag (pll N \. . . fil N *, a\ . . . af) + diag (xj, . . . , . . . 0) 

+ ^r;r. (5.2) 



Here I N denotes the N x N Identity matrix, x^ lies in the Lie algebra su(Nk) with 
x\ diagonal, M is given by 



K 



M = N-J^ N * 



and the tracelessness condition is 



M 



k=l 



K 



(5.3) 



(5.4) 



k=l 



k=l 



Note that the ordering of ps and as in f>.2\ is not significant; for example we could 
chose to take x% so that the elements are in decreasing order down the diagonal and 
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then degenerate and non-degenerate eigenvalues would be all mixed up in general. 
If x\'"' = then the eigenvalues of x± are exactly degenerate; we have K blocks of 
eigenvalues p\"' K with degeneracy Nk and singleton eigenvalues a\'"' M . When the 
exact degeneracy is relaxed slightly we have blocks of eigenvalues Aj^, m = 1, . . . N k 
with each block having central value 



k „k 



Ik 

m=l 



Pi (5-5) 



together with the singleton eigenvalues \^ + ^ = a{, j = 1, . . . M. 

For each sub-algebra H there is a contribution SU ^{R) to T^ 9 SU ^(R). This 
comes from the integration region V where the p\, a{, x\ are such that the eigenvalues 
of Xi satisfy 



V : 



\k _ \k 

m m' 



> e, k k', m = 1, . . . N k , m! — 1, . . . Ny 
\ k m -\ k m+1 \<e, m=l,...N k -l. (5.6) 



We note that any sequence of eigenvalues can be arranged in the manner implied 
by ^]6| for some H. Therefore by considering the Z|f su(N)(R) ^ or au possible H we 
exhaust all possible nearly degenerate eigenvalue configurations for x\. That is to 
say the degenerate term in Td,su(n)(R) is now the sum of contributions from all the 
possible Hs 

^D 9 SU{N)(R) = y^^D,SU(N) (R)- (5-7) 
H 

Now we bound the 1^ SU ^(R). The Vandermonde determinant for X\ is easily 
bounded from above through 

A= II ^m-^i') 2 II ( Xk m- Xk m>) 2 

k>k' ,m,m' k,m>m' 

<2 n * n 



k,m>m' 

K 



n H = N(N N *( N >° ~ !) ( 5 - 



k=l 



where we have used the fact that none of the eigenvalues can have magnitude more 



than 1 on account of the constraint |2.1| . Of course this bound is simply the product of 
the Vandermonde determinants for the constituent su(N k ) factors of the sub-algebra 
H . Now we note that the r", which are those off-diagonal elements of £2,...d that do 
not lie in H, are constrained by |2.1(J| and, following our procedure in the SU(3) case, 
we integrate them out to get 

r d / K m-i N 

zEsuww < c iR -^ d -^-v / n rK n h 

Jp »=1 \k=l j=l j 
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K / D \ / K M 

x ii /~n^ * h-E^x+^w-E^ 

fc=l V/** ^=1 / V k=l j=l 

(K M 
Tr 5*5* + JV fc pJpJ + £ ojoj - (l - n* (ei^) " 
fe=l i=i 

where cr^ is given by |5]4], the region lZ k is defined by 

[5*5*] <i?-( 2 -"), V/w,j. 



(5.9) 



(5.10) 



The right hand side of |5.9| is now bounded above by dropping the V constraint 
and integrating out the p k and which leaves us with 



K 



1-J> 



k=l 



(5.11) 



where 



fc=i 



D 



D,SU{N k ) 



(5.12) 



We note in passing that if H is empty then |5.11| simply reduces to the generic case 
2.12| . We now repeat the steps ^L7\ to fO| to find that 



JO 

The final step is by induction on iV: 
1- Fd,su{2) is given in fO|. 

2. From our results for SU(3) |4. 10 and 4.11 we deduce that 

Fd,su® < C 3 R- 4 ^ D , D > 4 



(5.13) 



(5.14) 



3. For gauge group S77(4) the possible sub-algebras, if, are sw(2), sw(2) © sw(2) 
and sw(3). For D = 3 we find that 



X 3 ,5C/(4)(i?) <C 5j R" 4S log i? 



(5.15) 



and hence the integral for Z$ t su(<i\ converges. For D > 4 it is simple to check 
that Zd,suU) converges and the dominant term in T d ^ S u(a){.R) comes from 
H = su(3) (we will give a general formula for the behaviour of 1d,su(n){R) 
below). 
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4. We now assume that convergence is established for D > 3 for iV < N*. As a 
consequence we have that 

Fd,su{N) < CR- (2 - ri)D{N2 ~ 1)/2 , D>3, and 3 < iV < N* (5.16) 

together with the slightly different bounds [4.9| and |5.14| for N = 2 and 3. This 



is enough information to bound Z|f SU ^ N ^(R) using |5.11| for any sub-algebra H 

lDSU{N*)(R) < C^'^-^+^^iNl-l)) [ R n 2 6 D , 3 ( hgR yn 2 6 D , 3+ n 3 6 D ^ 

(5.17) 

where n 2 ^ denotes the number of su(2) and su(3) factors respectively in H. It 
is straightforward to check that for D > 3 the slowest decay at large R occurs 
when H = su(N-l); basically this minimises the number of R~ 2 factors coming 
from off-diagonal elements not in the sub-algebra and maximises the number 
of i? -1 factors coming from elements in the sub-algebra. Thus for D > 3 and 
iV > 4 we have 

X d 9 su(n)( R ) < CR-W-W-V+wW-Q^logR) 6 "*** ' 3 (5.18) 



Applying the criterion [2.5| shows that both generic [2.12| and degenerate terms 
make a finite contribution to Zd,su(n) for D > 3, N > A. 

This completes our proof. 



6. Discussion 



6.1 Correlation Functions 



We can extend the definition of the partition function |1.1| to correlation functions so 
that 



< . > 



D poo ( \ 

J] / dX M (.)exp PTTr [X„X U ] 2 ) 

^ =1 J-oo \^ >u J 



(6.1) 



where (.) represents some kind of product of the with P factors. Making the 



change of variables |2.1| , realising that the absolute value of the corresponding product 
of the must be bounded by a constant, and using |5.18j we get that 



R 



<.>\<C, dRR- 2N(D -V +3D - 5 R p . 



(6.2) 



Thus correlators with fewer than 



P c = 2N{D - 2) - 3D + 4 



(6.3) 
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factors are guaranteed to be finite; of course correlators with more factors than this 
may be finite but then they must have some special property so that the leading di- 
vergences cancel. The authors of ]14j "guessed" on the basis of reasonable arguments 
that the eigenvalue density p(A) for behaves like 

P (A)~ Const \-mD-2)+3D-5 (6 4) 



at large A. This is completely consistent with our results. 



6.2 Supersymmetric Integrals 

The supersymmetric partition functions are given by 



Z d% = II / dX » V ^,g(X,) exp Tr ^ X ^ 



(6.5) 



where the Pfaffian Vd,g arises from integrating out the fermionic degrees of freedom 
and is a homogeneous polynomial of degree (D — 2)(iV 2 — 1). We can of course 
regard Z^f G as being a correlation function in the bosonic theory and apply the 
considerations of section 6.1 to it for P = (D — 2)(N 2 — 1); we find immediately 
that all the supersymmetric partition functions are naively divergent. However the 
Pfaffian contains many terms with different signs and we expect many cancellations. 
The simplest example is for SU(2) where the Pfaffians are known explicitly [12| and 
(except for D = 3 where Z^ s su( ^ = because T > ^ ) su(2) is an odd function) can be 
expressed as sums of powers of commutators [Xp, X v \. This is particularly convenient 
with our method because the rescaling in |2.1| followed by the restriction to the region 
1Z\ in p.6| means that we can bound Zffsu(2) simply by setting all commutators to 
a constant. It follows that if the bosonic partition function converges so does the 
supersymmetric one; thus the D = 6 and 10 partition functions are convergent but 
D = 4 is marginal and we would need to work harder (in fact it is known to converge). 
The situation with bigger N is more complicated mainly because relatively little is 
known about the Pfaffians and we will return to this problem in a separate paper. 



6.3 Other Gauge Groups 

In this paper we have concentrated on SU (N) gauge groups. However all the results 
we obtain depend in a well defined way on simple group theoretical properties such as 
the order, rank and sub-algebras of Lie algebras. It is therefore tempting to suppose 
that expressed in this form our results would carry over directly to any Lie group. 
However some of the steps we have made such as the diagonalization of X\ and the 
inductive argument in section 5 do depend on the group being SU(N) and we will 
deal with the other groups in a separate paper. 
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A. Divergent Matrix Integrals 

The D = 2 integral 

/oo 
dX lC /X 2 exp(Tr [X 1; X 2 ] 2 ) (A.l) 
-oo 

is divergent for all SU(N). This is easily seen by diagonalizing X\\ the integrand 
then does not depend upon the diagonal elements of X 2 and so the integral over 
them diverges. 

Some other low N and low D integrals are divergent. To see this we go back to 



TTT| , diagonalize Xi, and separate out the diagonal elements for v > 1, 

X u = di&g{\ ul ,...\ uN )+X^. (A.2) 

We then change variables from the X„ to (D — l)(g — l) dimensional polar coordinates 
with radial variable r and angular variables {Oi}. The integral over all the diagonal 
elements is gaussian so we do it and are left with 

/OO P 
r D( 9 -2i)-2 g+ 2i-i dr J ^ Fl ( {e .})exp(-r 4 F 2 ({^})) (A.3) 

where F\ and _F 2 are horrible but positive semi-definite functions and Q is the (D — 
l)(g — l) dimensional solid angle. We see immediately that the integral over r diverges 
at r = if 

D < ^ (A .4) 
g-% 

so we deduce that the D = 3, 4 integrals diverge for SU (2) and that the D = 3 
integral diverges for SU(3). 
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